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FUNCTIONALS OF DIRICHLET PROCESSES, THE 
CIFARELLI REGAZZINI IDENTITY AND BETA-GAMMA 

PROCESSES 1 

By Lancelot F. James 

Hong Kong University of Science and Technology 

Suppose that Pg(g) is a linear functional of a Dirichlet process 
with shape 6H , where > is the total mass and H is a fixed proba- 
bility measure. This paper describes how one can use the well-known 
Bayesian prior to posterior analysis of the Dirichlet process, and a 
posterior calculus for Gamma processes to ascertain properties of 
linear functionals of Dirichlet processes. In particular, in conjunction 
with a Gamma identity, we show easily that a generalized Cauchy- 
Stieltjes transform of a linear functional of a Dirichlet process is 
equivalent to the Laplace functional of a class of, what we define 
as, Beta- Gamma processes. This represents a generalization of an 
identity due to Cifarelli and Regazzini, which is also known as the 
Markov-Krein identity for mean functionals of Dirichlet processes. 
These results also provide new explanations and interpretations of 
results in the literature. The identities are analogues to quite useful 
identities for Beta and Gamma random variables. We give a result 
which can be used to ascertain specifications on H such that the 
Dirichlet functional is Beta distributed. This avoids the need for an 
inversion formula for these cases and points to the special nature of 
the Dirichlet process, and indeed the functional Beta-Gamma calcu- 
lus developed in this paper. 

1. Introduction. Let P denote a Dirichlet random probability measure 
on a Polish space y, with law denoted as T>(dP\9H), where 9 is a nonnega- 
tive scalar and H is a (fixed) probability measure on y. In addition, let A4 
denote the space of boundedly finite measures on y. This space contains the 
space of probability measures on y. The Dirichlet process was first made 
popular in Bayesian nonparametrics by Ferguson (1973) [see also Freedman 
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(1963) for an early treatment], and has subsequently been used in numerous 
statistical applications. Additionally, the Dirichlet process arises in a vari- 
ety of interesting contexts outside of statistics. Formally, P is said to be a 
Dirichlet process if and only if for each finite collection of disjoint measurable 
sets Ai, . . . ,Ak, the random vector P(A\), . . . ,P(Ak) has a Dirichlet distri- 
bution with parameters 0H(A\), . . . ,0H{Ak). In particular, P(A) is a Beta 
random variable for any measurable set A. An important representation of 
the Dirichlet, which is analogous to Lukacs characterization of the Gamma 
distribution, is P(-) = //(-)/T where fx is a Gamma process with finite shape 
parameter OH and T = fy fj,(dy) is a Gamma random variable with shape 
and scale 1. The law of the Gamma process is denoted as Q(d/J>\QH) and is 
characterized by its Laplace functional 



for each positive bounded measurable function g on y. For our purposes we 
shall consider the more general class of real- valued functions g which satisfy 
the constraint 



Jy 

This condition, (1), as shown by Doss and Sellke (1982) and Feigin and Tweedie 
(1989), is necessary and sufficient for the existence of the linear functionals 



An important fact is that T and P are independent, which as we shall see, 
has a variety of implications. 

An interesting problem initiated in a series of papers by Cifarelli and Regazzini 
(1990) is the study of the exact distribution of linear functionals P(g) of the 
Dirichlet process. One of their contributions is the important identity 



where typically z is in the complex plane c & . We call (2) the Cifarelli- 
Regazzini identity. The result in (2) is sometimes called the Markov-Krein 
identity for means of Dirichlet processes. Diaconis and Kemperman (1996) 
discuss some consequences of this result, which has implications relative to 
the Markov moment problem, continued fractions theory, exponential repre- 
sentations of analytic functions, and so on [see Kerov (1998) and Vershik, Yor and Tsilevich 
(2001)]. Vershik, Yor and Tsilevich (2001) expand upon this, emphasizing 
that the right-hand side of (2) is the Laplace functional of a Gamma pro- 
cess with shape OH. That is, 




(1) 




P(g) = I g{y)P{dy) 




(3) 
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Their interpretation, which is in the sense of a Markov-Krein transform, is 
that the generalized Cauchy-Stieltjes transform of order 9 of P(g), where P 
is a Dirichlet process with shape OH, is the Laplace transform of fi(g) when 
fi is the Gamma process with shape 9H. The authors then exploit this fact 
to rederive (3) via an elementary proof using the independence property of 
P and T. An interesting question is, what can one say about 

when 9 and q are arbitrary positive numbers? That is, can one establish 
a relationship of (4) to the Laplace functional of some random measure, 
say /j,*, which is similar to //, for all q and 91 Lijoi and Regazzini (2004) 
establish analytic results for (4), relating them to the Lauricella theory of 
multiple hypergeometric functions. Theorem 5.2 of their work gives ana- 
logues of (2), stating what they call a Lauricella identity, but does not 
specifically state a relationship such as (3). We should say for the case 9 > q 
that it would not be terribly difficult to deduce an analogue of (4) from 
their result. However, this is not the case when 9 < q, which is expressed in 
terms of contour integrals. Their representations, for all 9 and q, as clearly 
demonstrated by the authors, however have practical utility in regards to 
formulae for the density of P(g). In this case, one wants to have an expres- 
sion for (4), when q = 1 and for all 9. Some related works include the papers 
of Kerov and Tsilevich (1998), Regazzini, Guglielmi and Di Nunno (2002), 
Regazzini, Lijoi and Priinster (2003) and the manuscript of James (2002). 

1.1. Preliminaries and outline. In this paper we develop results that are 
complementary to the work of Lijoi and Regazzini (2004) and Vershik, Yor and Tsilevich 
(2001). In particular, we show that (3), as interpreted in Vershik, Yor and Tsilevich 
(2001), extends to a relationship between (4) and the Laplace functional of a 
class of what we call Beta-Gamma processes defined by scaling the Gamma 
process law by T~ d , for all numbers d such that 9 — d > 0, that is, processes 
with laws equal to 

(5) BG(d^\9H, d) = J^^T- d Q{dii\0H). 

In particular, our main result concerns the choice of d = 9 — q for arbitrary 
positive numbers which are not necessarily equal. The approach relies in part 
on, in this case, mostly familiar Bayesian prior posterior calculus for Dirichlet 
and Gamma processes in conjunction with the usage of the following well- 
known Gamma identity for q > 0: 



(6) 
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That is to say, purely analytic arguments are replaced by Bayesian argu- 
ments using the familiar results in Ferguson (1973), Lo (1984) and Antoniak 
(1974), thus giving the derivations a much more interpretable Bayesian fla- 
vor. More specifically, albeit less well known, we use the results in Lo and Weng 
(1989) as demonstrated for more general processes in James (2002). This 
bypasses the need, for instance, to verify certain integrability conditions 
and the usage of limiting arguments. Moreover, somewhat conversely to 
Lijoi and Regazzini (2004), we show how properties of the Dirichlet and 
Beta-Gamma processes yield easily interesting identities related to Lauricella 
and Liouville integrals [see Lijoi and Regazzini (2004) and Gupta and Richards 

(2001) ]. Although we exploit the independence property of T and P to 
prove our results, our approach is quite different from the methods used by 
Vershik, Yor and Tsilevich (2001) to prove (3). While their proof is certainly 
elegant, it does not seem possible to extend to other processes. Our methods, 
however, are influenced by their proof of an analogous result for the two- 
parameter extension of the Dirichlet process [see Pitman (1996)] which relies 
on (6) and the fact that such processes are based on scaled laws. That is to 
say, we present an approach which is extendable to other models [see James 

(2002) , Section 6]. However, for the Beta-Gamma processes defined in (5), 
the independence property between T and P translates into the property 

(7) / h(P)V(dP\9H)= [ h(P)g(dfi\6H)= [ h{P)BG(dfi\9H,d) 
Jm Jm Jm 

for all integrable h. The property (7) seems to suggest that the Beta-Gamma 
process may not have much utility relative to calculations involving P; how- 
ever, it is precisely this property that we shall exploit. In the next sec- 
tion we shall first develop, rather quickly, two supporting results concerning 
the calculus of Dirichlet and Beta-Gamma processes. We will then show 
how these results are used to easily derive our main results in Theorems 
2.1 and 2.2 based on Bayesian arguments. We close the paper by showing 
how our methodology, a Beta-Gamma calculus for Dirichlet processes, leads 
to a functional analogue of the classical Beta-Gamma calculus for random 
variables. That is, we provide conditions on H such that Pe(g) is Beta dis- 
tributed. 

2. Functional of Dirichlet processes, the Cifarelli Regazzini identity and 
Beta-Gamma processes. We start by recalling some properties of the Dirich- 
let process. Let Yi,...,Y n denote random elements in the space y, which 
conditional on P are i.i.d. with law P. P is a Dirichlet process with shape 
9H. These specifications define a joint law of (Y,P), where, from Ferguson 
(1973), it follows that the posterior distribution of P\Y is also Dirichlet with 
shape 

n 

(9 + n)H n = eH + J2^Y l - 
i=i 
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Additionally, the marginal distribution of Y is the Blackwell and MacQueen 
(1973) distribution described as 

ndY\9H) = jy^OHidYi) f[ (eH + £ Sy^j (dY). 

The Blackwell-MacQueen distribution admits ties. Hence one can represent 
Y = (Y*,p), where Y* = {Y{ , . . . , Y*^} denotes the n(p) < n unique values 
of Y. The expression p = {C\, . . . , C n ( p )} denotes a partition of the integers 
{1, . . . , n} with n(p) elements. The Cj = {i:Yi = Y*} for j = 1, . . . , n(p) 
denote the collection of values equal to each unique Y?, for j = 1, . . . , n(p). 
The cardinality of each cell Cj is denoted as ej n . When H is nonatomic, 
then one can write 

n(p) 

F{dY\6H) = n(p\0) ]J H(dY*), 
i=i 

where 

is a variant of Ewens' (1972) sampling formula, which was independently 
derived by Antoniak (1974). It is also called the Chinese Restaurant pro- 
cess [see Pitman (1996)] and plays a fundamental role in Lo (1984). If H is 
discrete with probability mass function p, then 

T(0) n ^T(e P (Y*) + e jn ) 

In any case, note that appealing to standard Bayesian arguments, the results 
of Ferguson (1973) imply that one has 

(8) / h{P)V{dP\6H)= ( [ h(P)V(dP\(9 + n)H n ) F(dY\9H). 
Jm Jy n Um J 

This simple consequence is fundamental to our presentation. It is evident 
that (8) along with the various forms of P(dY\6H) yield nontrivial expres- 
sions which might otherwise require an appeal to, for instance, the theory 
of special functions or combinatorics. In the same spirit, we now derive the 
posterior distribution of the Beta-Gamma processes. From Lo and Weng 
(1989), one has the following disintegration of measures: 

(9) n{dYi)G{dii\eH) = r( p ( + n) g \dfi\9H + J2 S Y ]ndY\9H), 
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where Q(dfj,\6H + Ya=i ^vj denotes a Gamma process with shape 9H + 
Er=i% Usin g ( 7 )> ^ is easy to see that J M U?=i P(dYi)Bg{dfi\9H, d) = 
¥*{dY\9H). Combining this fact with (9) easily yields the following descrip- 
tion of the posterior distribution of a Beta-Gamma process. 

PROPOSITION 2.1. Let fi have law Bg(dfi\9H,d) defined for all d, such 
that 9-d>0. Then from (7), the law of P = n/T is V(dP\9H). Suppose 
that Yi,...,Y n \P are i.i.d. P; then the posterior distribution of fi\Y is a 
Beta-Gamma process with parameters (9 + n)H n and n + d, defined as 

Bg(dfj,\(e + n)H n , n + d) = ^ + n J T-^g(d^\(9 + n)H n ). 

Hence, similar to (8), one has 

f{n)Bg{dn\6H,d) 



(10) 



M 



f(fj,)Bg(dfi\(9 + n)H n ,n + d) 

M 



F(dY\9H) 



for all integrable f . Note that setting d = shows that if \i is g{d[i\9H) , 
then its posterior distribution is Bg(dfi\(9 + n)H n ,n) , which is not a Gamma 
process. 

Remark 2.1. Note that the use of (7), (8) and (10) sets up a myriad of 
interesting equivalences which will prove useful in our derivations. However, 
we do point out that while (10) implies (8), the converse is not true. 

Another important property of the Gamma process that we shall exploit 
is the algebraic identity 

(11) / e-^ T+w ^))g {MeH) = {l + v) -e e -J y ^[i+M(i+v)) 9 (y)]eH(dy)_ 
JM 

Let 

B(du\a, b) = £'/ ^//, « a " 1 (l " «) 6_1 du for < n < 1 
T(a)T(b) 

denote the density of a Beta random variable with parameters (a,b). We 
now establish our final preliminary result before our main theorem. 

Proposition 2.2. Let 9 and q be arbitrary nonnegative numbers. Then 
for any integer n > that satisfies the constraint 9 + n — q > 0, the following 
formula holds: 

T(9 + n) f 1 



/ ^S{dfi\(9 + n)H n 
T(9 + n- q) J M (T + zMfiO) 

(12) , 

= / e~ z ^Bg(dfi\{9 + n)H n ,9 + n - q) 
JM 
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Proof. Apply the Gamma identity to (T + /j>(g))~ q and then (11) with 
v = w and 9H replaced by (9 + n)H n , to show that the left-hand side of (12) 
is equal to 

r ( g + n ) f°° v q-\~ f y lo S [l+(v/(l+v))zg(y)]6H(dy) 



T(9 + n-q)T(q) 



n(p) . \ — e- 



Similarly, the following expression is obtained for the right-hand side of (12) 
by applying the Gamma identity to x~^ 9+n ~ q \ (11) with w = 1 and a further 
change of variable: 

(13) / e -fy^+uz 9 (y)]0H(dy) jt {1 + uzg( Y*))-^J3(du\q,e + n-q). 
Jo 3=1 

The result is obtained by applying the transformation u = v/(l +v). □ 

We now present a new result which relates the generalized Cauchy-Stieltjes 
transform of Dirichlet process linear functionals to the Laplace functional 
of Beta-Gamma processes. This presents a generalization of the Cifarelli- 
Regazzini identity, complementary to the Lauricella identities deduced in 
Lijoi and Regazzini [(2004), Theorem 5.2]. We also present some interesting 
additional identities. 

Theorem 2.1. Let T)(dP\6H) denote a Dirichlet process with shape 9H . 
Let g denote a function satisfying (1). Then the following relationships are 
established: 

(i) For any positive q and 9, 

(14) / {l + zP{g))~ q V{dP\9H)= f e- z ^ 9) Bg(dfi\9H,9-q). 
Jm Jm 

Note that the law BG(dfi\9H, 9 — q) exists for all positive 9 and q, and arbi- 
trary H, since 9 — (9 — q) = q>0. 

(ii) For any positive q and 9, and any integer n > which satisfies 9 + 
n — q > 0, the quantities in (14) are equivalent to fy n [f M e~^ 9 ^ BG(dfx\(9 + 
n)H n ,9 + n — q)]P(dY\9H). An explicit expression can be deduced from the 
equivalence of the inner term to (13). In particular, when H is nonatomic, 
the expression is equivalent to 

Y 7T(p|0) f 1 e'Sy ^ll+uzg(y)]eH(d y ) 

„ Jo 



™(p) „ 

n / (l + uzg(y)y e ^H(dy) 

3=1 Jy 



B(du\q, 9 + n — q). 
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For the Gamma process with law Q(dfj,\9H), its Laplace functional may be 
represented as above for all n > 1 and q = 9. 

(iii) When 9 — q > 0, statements (i) and (ii) with n = imply that 

( (l + zP{g))- q V(dP\9H) 
JM 

( 15 ) .1 f 

e -Jy^[l+uz 9 (y)]eH(dy) B{du ^ e _ q ^ 





which coincides with the result in Lijoi and Regazzini [(2004), Theorem 5, 
equation (5.2)]. 

Proof. A general strategy is formed by first writing (1 + zP{g))~ q = 
T q (T + zfi(g))~ q = h(P). For the proof of statement (i), we first assume 
without loss of generality that q = n + d, where d is a positive number such 
that 9 — d > 0, and n > is an integer chosen such that 9 + n — q > 0. This 
means that T« = T n+d . Now using (8) and then (7) with BG(dfi\(9 + n)H n , q) 
yields 



I {l + zP{g)r q V{dP\9H) 
JM 



IM 

T(9 + n) 



d)Jy 



M iT + z^g)) q 



G(dfi\(9 + n)H n ) 



¥(dY\9H). 



Apply Proposition 2.2 to the inner term, recalling that 9 + n — q = 9 — d. 
This yields the desired expression, 



(16) 



/ [/ e- z ^Bg{d^\{9 + n)H n ,9 + n 
= [ e- z ^BQ{d[i\9H,9-q). 

JM 



F(dY\9H) 



Note how again an appeal to a Bayesian argument, that is, using (10) in 
Proposition 2.1 with /(//) = e~ z ^ 9 \ is used to deduce easily the equivalence 
of the right- and left-hand sides of (16). □ 

Remark 2.2. Since H is an arbitrary distribution, the result applies 
to a Dirichlet process posterior distribution based on, say, a sample of size 
m having no particular relationship to n. For concreteness, suppose P is 
a Dirichlet process with shape aPo + Y^Li > where Po is an arbitrary 
probability measure, a is a positive scalar and X\, . . . ,X m are fixed points. 
Then the results in Theorem 2.1 hold for this P by setting 9 = a + m and 
9H = aP + YZi8 Xi - 
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Remark 2.3. As in Kerov and Tsilevich (1998) and Vershik, Yor and Tsilevich 
(2001), Theorem 2.1 applies to the joint distribution of linear functionals, 
say (P(gi),...,P(g m )), where g 1 ,...,g m are functions satisfying (1). The 
generalized Cauchy-Stieltjes transform for joint distributions is defined by 



replacing zP(g) by YT=i *iP(9i)- Since YT=i *iP(9i) = P{YT=i *i9i), the re- 



suit is seen by replacing zg with Y^=i z i9i i n Theorem 2.1. 

We now discuss some interesting results obtained from Theorem 2.1. Note 
the relative ease by which Bayesian arguments can be used to derive oth- 
erwise complex expressions such as that appearing in Theorem 2.1(h). The 
case q = 1 is of particular interest in terms of giving an expression for the 
Cauchy-Stieltjes transform P(g), which can be inverted to obtain an ex- 
pression for the distribution of P(g). Setting q = 1 in Theorem 2.1(h) gives 
a variety of equivalent expressions which hold for all 9 and n > 1. Here, as a 
corollary, we present the simplest expression that holds for all 9 with n = 1. 

Corollary 2.1. Let D(dP\9H) denote a Dirichlet process with shape 
9H , where H is an arbitrary probability measure. Let g denote a function 
satisfying (1); then for all 9 > 0, 



When 9 > 1, this expression equates to the expression in Theorem 2.1 (hi). 
When 9 = 1, Theorem 2.1(h) shows that the right-hand side of (17) is the 
Laplace functional of a Gamma process with shape H, which corresponds 



The expression (17) can be seen as complementary to the expressions ob- 
tained in Lijoi and Regazzini (2004). However, our results are quite different 
in the case where < 9 < 1, where those authors obtained an expression in 
terms of contour integrals. 

Let C(Z) denote the law of a random element Z. For the remainder of 
this work, let fJ,$ e—q be a Beta-Gamma process with parameters (9H, 9 — q), 
and let U a ,b denote a Beta(a, b) random variable. Let T a denote a Gamma 
random variable with shape a and scale 1, and let Y\ be a random element 
with distribution H . Let fig denote a Gamma process with shape 9H and 
assume that the variables fxg, U at b,T a ,Yi are independent. Additionally, let 
Pg denote a Dirichlet process with shape 9. When convenient we will simply 
write X = T to denote that the distribution of X is equivalent to that of T. 
That is, X = U a ^ means that X has a Beta distribution with parameters 




(17) 




1-u) 



to (2). 
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(a, b). The next result involves a series of distributional identities. These 
are based on Bayesian mixture representations deduced from the form of 
the posterior distribution mixed over the marginal distribution, P(dY\0H). 
Some important consequences will be demonstrated thereafter. 

Theorem 2.2. Let f-tg t g- q be a Beta- Gamma process with parameters 
(9H,9 — q) and let \ig denote a Gamma process with shape 9H. Then for 
all positive 9 and q and an integer n chosen such that 9 + n — q > 0, the 
following distributional equalities hold: 

(i) For all 9 > and q, and an integer n chosen such that 9 + n — q > 0, 

/ "(p) \ 

(18) C{^e,e-q) = £ I U q7 g +n - q fJ,g + U q fi +n _ q Gj tn 5y* J , 

where conditional on p the distinct variables on the right-hand side are mu- 
tually independent such that U q> g +n - q is Beta with parameters (q,9 + n — q), 
fig is a Gamma process with shape 9H , {Gj n } are independent Gamma ran- 
dom variables with shape e Jjn and scale 1. The distribution o/Y= (Y*,p) 
is W(dY\8H). In particular, if H is nonatomic, the Y* for j = 1, . . . , n(p) 
are i.i.d. H, and the distribution of p is 7r(p|0). Statement (i) implies the 
following results. 

(ii) For all 8 and q = 1, 

(19) C(Mo,e-i) = L(JJ\m + U lfi Ti8 Yl ). 

If fig denotes a Gamma process with arbitrary shape parameter 9H , then 

(20) C(flg) = C(U e ^e + Ug^SyJ. 
(hi) For all positive 9 and q, 

(21) C(fig^ q )=C(T g P e ), 

where T q is a Gamma random variable with shape q and scale 1 independent 
of Pq, which is a Dirichlet process with shape 9H . Hence for all positive 9 
and q, C{Tg^g^ q ) = C(T q fig), where Tq is Gamma with shape 9 and scale 1, 
independent of fJ,g t g- q . Similarly, T q and fig are independent. 

Proof. The distributional identity in (i) is a direct consequence of the 
mixture representation of the law of fig.0- q , in the form of the posterior 
distribution of fig g- q \Y and P(dY|0ff), deduced from the expression for the 
Laplace functional in Theorem 2.1(h). Note that all quantities on the right- 
hand side of (18), including p, are random. We now show statement (hi) 

follows from statement (i). Notice Tg +n : = ^g(y) + J2]=i Gj,n is a Gamma 
random variable with shape 9 + n independent of U q) g +n _ q . Moreover, using 



BETA-GAMMA PROCESSES 



11 



the mixture representation of the Dirichlet process derived from its posterior 
distribution and P(dY\9H), it follows that (fig + Y^=i Gj >n SY*)/Tg +n is a 
Dirichlet process with shape 9H, independent of Tg +n and U q fi +n _ q . Hence 
the right-hand side of (18) can be written as U q: g +n _ q Tg +n Pg. The result is 
completed by noting that U qi g +n - q Tg +n is equal in distribution to T q . □ 

Remark 2.4. The distributional identities in (19) and (20), which are 
new, are analogous to similar identities for Dirichlet processes which have a 
variety of applications, as can be seen in Diaconis and Kemperman (1996), 
Sethuraman (1994) and Hjort (2003). In addition, distributional results be- 
tween Beta and Gamma random variables have quite striking consequences, 
as can be seen from, for instance, Dufresne (1998). We view our results as 
functional extensions of some of those ideas, and it seems worthwhile to pur- 
sue more analogous results. Note importantly that our results do not require 
that fig(g) is Gamma distributed. 

Remark 2.5. The expression in Corollary 2.1 is obtained by evaluat- 
ing the Laplace transform of the right-hand side of (19), in the order of 
integration of fig, T%, Y\ and finally U\q. It is evident that other equiva- 
lent expressions can be formed by changing the order of integration. It is 
no coincidence that U\ g has the same distribution as T\/{T\ + T) where 
T = Jy fi{dy) = Tq. Additionally, further representations can be obtained by 
using the distributional identity 



where T\,Tg and t p are all independent and t p is a stable random variable 
with index < p < 1 . 

2.1. Distributional characterizations via the Beta- Gamma calculus. The 
expression (21) tells us precisely that, for all 9 and q, a Beta-Gamma process 
with parameters 9H and 9 — q is equivalent in distribution to a Dirichlet 
process with shape 9H , scaled by an independent Gamma random variable 
with shape q. Hence, using this interpretation the first result in Theorem 2.1 
is an immediate consequence of 



JM 

Although this viewpoint at first may seem to have limited usage, it has 
interesting consequences when combined with our other results, within the 



Ui, 
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context of the Beta-Gamma calculus. Note that since fj,g t g- q (g) = T q Po(g), 
one may apply the special features of the classical Beta-Gamma calculus 
combined with our results to deduce the following characterization of when 
Pg(g) is Beta distributed, in the case where g(Y) is not an indicator variable. 

Proposition 2.3. Let 0, q and a denote positive real numbers. Suppose 
that for some q > 0, /ig t g- q (g) = T a , that is, it is Gamma distributed with 
shape parameter a and scale \; then for the case < a < q, Pg(g) = U a ^ q - a 
and He(g) = TgU a:q - a . From the distributional identity (18), this is true if 
one chooses the distribution of g(Y) such that, for a fixed n > 1, 

n(p) 

(22) T e U a ^ a + J2 G 3 g(Y*) = T 6+n U a , q - a) 

i=i 

where TgU aA - a and J2j=i Gjg(Y*) are independent. For clarity, when n = l, 
(22) specializes to TgU aA -. a + T 1 g(Yi) = T e+1 U a)q - a . 

Proposition 2.3 provides a characterization for the reverse question as to 
which choice of H produces a Beta distribution for Pg(g). This is of course 
seen to be equivalent to specifying H to induce a particular distribution 
on the quantity Y^j=i GjgiY^), for some fixed value of n, which satisfies 
the constraint (22). It is clear by using (18), that this particular feature, of 

inducing a distribution on Y^j=i Gjg(Y J *) satisfying appropriate constraints, 
can be applied to any choice of distribution for Pe(g)- However, because of 
the available independence properties between Beta and Gamma random 
variables, the occurrence of a Beta distribution for Pe(g) can be checked 
several ways, not available to other distributions. In particular, note that 
Peid) = U a ^- a if and only if [iefi~ q = T a . Hence one can choose g(Y) (or 
check this), such that the Laplace transform 

E[ e -*i»' 8 -*to)] = (i + z )- a j 

or perhaps more easily using the Gamma process to check whether ne(g) 
satisfies 

e -f y log[l+zg(y)]eH(dy) = E r-zT e U a>g ^ a ] 

We show in the next proposition the limitations, within the context of Propo- 
sition 2.3, of choosing Y^=i Gjg(Yj) to be a Gamma random variable. 

Proposition 2.4. Suppose that H is chosen such that for some fixed 
n > 1, J2j=i Gjg(Y J *) has a Gamma distribution with shape parameter c n , 
depending on n, and scale 1. Then the only value of c n such that Pe(g) has a 
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Beta distribution is c n = n/2. Moreover, the distribution of Pe(g) is a sym- 
metric Beta distribution with parameters (6 + n/2,9 + n/2), for alln>\ and 
all 6 > 0. Equivalently £i0 ,-(0+n)(sO = T e+n/2 and fi e (g) = T e U e+n/2) e+n/2, 
and hence is never Gamma distributed. These specifications correspond to 
the choice of a = 6 + n/2 and q = 26 + n in Proposition 2.3. 

Proof. The proof is obtained by applying Theorem 2 of Dufresne (1998), 
combined with the constraints deduced from (22) in Proposition 2.3. It fol- 
lows that the only solution is given by 

TeUg +n / 2 ,8+n/2 + ^n/2 = ^0+n^9+n/2,0+n/2- ^ 

It is already known [see Cifarelli and Melilli (2000)] that if g(Y) = ^1/2,1/2 1 
the arcsine law, then P$(g) = t r e+i/2,6»+i/2- That is the case of n = 1 in Propo- 
sition 2.4. The case for n = 2 corresponds to 

T 6 U e+1>e+1 + [ P T{ + T 1 ]g(Y 1 ) + (1 - p)T{g(Y 2 ) = T e+2 U e+1>e+1 , 

where T[, T\ are independent exponential (1) random variables. Y\ and Y 2 
both have distribution H, but may be tied, p is a Bernoulli random variable 
with success probability 1/(^ + 1), corresponding to the case where Y\ = Y 2 , 
from the Blackwell-MacQueen urn scheme. It is not immediately clear how 
to choose H such that \pT{ + T 1 ]g(Yi) + (1 - p)T{g(Y 2 ) = T\. 
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